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Abstract We amend a recent dispersive analysis of 
the anomalous ij decay process ij —>■ 7r+7r“7 by the ef¬ 
fects of the 02 tensor meson, the lowest-lying resonance 
that can contribute in the Trty system. While the net 
effects on the measured decay spectrum are small, they 
may be more pronounced for the analogous ij' decay. 
There are nonnegligible consequences for the rj transi¬ 
tion form factor, which is an important quantity for 
the hadronic light-by-light scattering contribution to 
the muon’s anomalous magnetic moment. We predict 
total and differential cross sections, as well as a marked 
forward-backward asymmetry, for the crossed process 
77r“ —> TT~r] that could be measured in Primakoff reac¬ 
tions in the future. 
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1 Introduction 

The decay rj —^ 7r+7r“7 is one of the processes driven by 
the chiral anomaly mm- The reduced scalar decay am¬ 
plitude (to be defined below) in the SU(3) chiral limit 
and at vanishing momenta is given entirely in terms 
of the electromagnetic coupling e and the pion decay 
constant 

^ aVStt^F^ ^ 5.65 GeV (1) 

Higher-order corrections to the anomaly can be eval¬ 
uated in chiral perturbation theory [3], and pion-pion 
rescattering in the final state resummed effectively us¬ 
ing dispersion theory [1]. Besides thus being an inter¬ 
esting decay in its own right to test our understand¬ 
ing of the interaction of light pseudoscalar mesons with 


photons, this decay is particularly noteworthy as a fun¬ 
damental ingredient in a dispersive analysis of the rj 
transition form factor rj — >■ 77* [^. This quantity is a 
crucial input necessary for the ongoing program to ana¬ 
lyze the hadronic light-by-light scattering contribution 
to the anomalous magnetic moment of the muon, com¬ 
bining as many pieces of experimental information as 
possible in a model-independent fashion laiz] . A simi¬ 
lar analysis has also been pursued for the transition 
form factor [8]. 

As pointed out in Ref. [1], the decays of 77 and rj' 
into 7 r“'" 7 r “7 pose a beautiful and simple example to 
demonstrate the universality of final-state interactions. 
Neglecting (tiny) contributions of F and higher partial 
waves for the pion pair, the authors show that the re¬ 
duced decay amplitude can be written as 

F{s,t,u) = P{t)F^{t), (2) 

where t = is the squared invariant mass of the pion 
pair, F^ (t) denotes the pion vector form factor as mea¬ 
sured in e^e“ —>• tt+tt", and P(t) is a polynomial. Com¬ 
parison to experimental data obtained by the WASA- 
at-COSY [5] and KLOE [TU] Collaborations demon¬ 
strated that within experimental accuracy, the poly¬ 
nomial can be assumed to be linear, P{t) = A{\ + at), 
with [TU] 

a= (1.32±0.13)GeV^. (3) 

This result gives rise to several interesting questions. 
Obviously, Eq. is only an approximation, tested suc¬ 
cessfully in the physical decay region, AM^ <t< M^. 
The universality of final-state interactions expressed 
therein is only valid in the region of elastic pion-pion 
rescattering, which is phenomenologically a good ap¬ 
proximation up to roughly t ft! 1 GeV^. From generic 
considerations about the asymptotic behavior of the 
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decay amplitude, one would rather expect P{t) to be¬ 
come constant for large t, such that the decay ampli¬ 
tude falls like 1 jt similar to the asymptotic behavior of 
{t). The continuation beyond the physical regime is 
interesting in particular with regard to the application 
within a dispersive integral to obtain the rj transition 
form factor [^, as in principle that integral covers all 
energies. 

The present article is built on the following observa¬ 
tion. If we continue the amplitude ([2]) naively to nega¬ 
tive t, we ought to observe a zero at or near t = — 1/a Ri 
—0.76 GeV^. Such a kinematical regime is indeed acces¬ 
sible: in the crossed reaction 77r“ —>■ TT~r], which could 
be measured in a Primakoff-type reaction, i.e., the scat¬ 
tering of a charged pion in the strong Coulomb field of 
a heavy nucleus, producing an additional ij. Such a Pri- 
makoff program is currently pursued by the COMPASS 
Collaboration (see, e.g.. Ref. [IT] for an overview), using 
a 190 GeV tt~ beam and cutting on very small momen¬ 
tum transfers in order to isolate the photon-exchange 
mechanism from diffractive background. In this way, 
COMPASS can investigate reactions to various fi¬ 
nal states, in particular Compton scattering in order to 
extract the charged-pion polarizabilities [HKTS], 
to investigate the chiral anomaly [T4l[T5], or three pi- 
ons testing chiral predictions [reifTT]. In this paper, we 
want to provide the theoretical motivation to also mea¬ 
sure the final state TT~r]^ as well as a prediction for the 
cross sections that are to be expected. 

For this purpose, beyond using crossing symmetry, 
we need to amend the amplitude m for the follow¬ 
ing reason. The assumption underlying Eq. © is the 
neglect of so-called left-hand cuts: the two pions under¬ 
going final-state interactions are assumed to originate 
from a point source, such that the amplitude is of form 
factor type, and any interaction (resonant or nonreso¬ 
nant) in the irr] channel is neglected. This approxima¬ 
tion can be justified at low energies by appealing to 
chiral perturbation theory: the Trrj P-wave is chirally 
suppressed (as well as all higher partial waves) [T51[T1] , 
an imaginary part only appears at three-loop order, any 
phase shift is therefore expected to be very small. Fur¬ 
thermore, the irr] P-wave has exotic quantum numbers 
= 1 and the search for possible resonances in 
this channel is not fully conclusive so far The 

first well-established resonance that is therefore going 
to be important in the process 'yir —^ irr] is the P-wave 
tensor meson 02(1320). To investigate its influence is 
important for several reasons: 

— its inclusion will demonstrate to what extent the 
feature expected from Eq. a zero (or at least a 
pronounced minimum) in certain differential cross 


sections, can survive in a more complete description 
of the amplitude; 

— it will provide a characteristic breakdown scale in 

the mj invariant mass squared s = above 

which Trrj resonances dominate the cross section; 

— finally, we can use the 02 as the likely most im¬ 
portant left-hand-cut structure for the decay 77 —?► 
7r+7r“7, to study to what extent it affects the decay 
amplitude, and whether its effect is consistent with 
the experimental decay data available. 

The outline of this article is as follows. In Sect. [2j 
we recapitulate the dispersive representation of the 
■q —>■ 7r+7r“7 decay amplitude of Ref. [3], before cal¬ 
culating contributions of the 02 tensor meson first at 
tree level, then including pion-pion rescattering effects 
dispersively. Section |3] compares the resulting observ¬ 
ables to the measured q —> Tr+Tr^y decay spectrum and 
briefly discusses the possible impact on the q transition 
form factor. In Sect. |4l we give our predictions for the 
crossed process 77r“ —>■ '!T~q, discussing total and differ¬ 
ential cross sections, the leading partial waves, as well 
as the resulting pronounced forward-backward asym¬ 
metry. We close with a summary. A brief discussion of 
the related decay q' —>■ 7r+7r“7 is relegated to an ap¬ 
pendix. 


2 ?7 —>■ 7r7r7 with left-hand cuts 

2.1 Amplitude, kinematics 


We write the decay amplitude for the process 

q{q) ^ Tr+{pi)Tr~{p 2 )j{k) (4) 


in terms of a scalar function F{s, t, u) according to 

M{s, t, u) = ie^vapt^{k)p'{p 2 q^t, u), (5) 

with the Mandelstam variables given as s = {q — Pi)^, 
t = {Pi + 7 ' 2 )^j and u = {q — p2)^. F{s, t, u) in the chiral 
limit fulfills the low-energy theorem 7^(0, 0,0) = F/ 7 r 7 i- 7 - 
The cosine of the t-channel center-of-mass angle is given 
by 


cos 6t = zt = 


s — u 





( 6 ) 


The t-channel partial-wave expansion is of the form 


J^(s, t, u) — Y. a'(=!.)A(!). (7) 

odd I 


where Pl{zt) denote the derivatives of the standard Leg¬ 
endre polynomials. Due to the strong suppression of F 
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Fig. 1 Graphical illustration of the discontinuity equa¬ 
tion (ED- The gray circle denotes the i-channel P-wave pro¬ 
jection of the Tj 7r+7r~7 decay amplitude, whereas the white 
circle stands for the P-wave pion-pion scattering amplitude. 


and higher partial waves at low energies, we will al¬ 
most exclusively be concerned with the P-wave, which 
is obtained by angular projection according to 


The differential decay rate with respect to the pion- 
pion invariant mass squared is given by 


d/" ^ . X 3 

^ = 4 


dzt(l - 


= Po(t)x (|/i(t)| 2 + ...), 

Po(i)= ^ 


12(87rM^)3 ’ 


(9) 


where the ellipsis in the second line represents neglected 
higher partial waves. 

In the absence of left-hand cuts and ignoring inelas¬ 
ticities, the P-wave should obey the following represen¬ 
tation [?]: 


fi{t) = P{t)n{t), f2(t)=exp|- f dx 

I JiM^ -t) } 

( 10 ) 

where Q{t) is the Omnes function [22] given in terms of 
the pion-pion P-wave phase shift 5{t) = and P{t) 
is a polynomial. The representation ca is a solution 
to the discontinuity relation 

disc/i(t) = (11) 

as obtained from elastic pion-pion rescattering; see 
Fig. [T] It obviously fulfills Watson’s final-state inter¬ 
action theorem m- the phase of fi{t) agrees with the 
elastic scattering phase S{t). In the following, we will 
take 6{t) from the representation given in Ref. [24] . 
As already pointed out in the introduction, comparison 
with data [siiin] suggested that the polynomial P{t) is 
linear in the decay region. 



Fig. 2 Tree-level contributions of the a2(1320) resonance to 
r] 7r+7r“7 in the s- (left) and u-channel (right). 


to very good accuracy. In fact, in Ref. |3], the Omnes 
function was replaced by the pion vector form factor 
(t), which is a phenomenologically attractive repre¬ 
sentation insofar as the latter is itself directly experi¬ 
mentally observable. Both representations are equiva¬ 
lent modulo a moderate shift in the parameter a ^ an 
due to the observation that the form factor is in turn 
proportional to the Omnes function up to a linear 
polynomial below 1 GeV, with a slope of the order of 

0.1 GeV^ [5]. 


2.2 Tree-level contribution of the 02(1320) 


We begin by calculating the tree-level contribution of 
the 02 tensor meson to the amplitude 77 —>■ 7 r“'" 7 r “7 as 
shown in Fig. [21 For the formalism of coupling tensor 
mesons to Goldstone bosons, we follow Ref. [55]. The 
single necessary interaction term required to describe 
the decay of a tensor meson into two pseudoscalars is 
given by 


Ctpp = gT{Tfi„{u^, 


(13) 


where (.) denotes the trace in flavor space. For simplic¬ 
ity, we only display the nonstrange SU(2) part of the 
tensor field relevant to our calculation explicitly, 


. _ J_ / al \/2a} \ 

V2[V2a^-a°2j 


(14) 


The Goldstone bosons are encoded in the field = 
i{u^dfj,u — ud^u^) (neglecting external currents), where 


( * 

u = exp I 


V2F^ 


•\/27r“ —TT^ -\- — 


V3 , 


(15) 


From Eq. m, we can calculate the decay width for 
02 —>■ 7177, employing the polarization sum of the 02 [55] 

£711^(0epcr(0 = 

pol 

Ppu,piy(}) = -^{PppPua + PvpPprr) ~ -^PpvPpa, 

p - a - hijL 

™a2 


P{t) = A(1 -I- ant) 


( 12 ) 


(16) 
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and find 
r{a2 -)> TTT]) 


ISOttF^ 


(17) 


where A(a, 6, c) = — 2{ab + ac + be) denotes 

the usual Kallen function. Equation (II3, with the total 
width Eaj = (107 ± 5) MeV and the branching fraction 
B{a 2 —>■ Trrj) = (14.5 ± 1.2)% [55], leads to the coupling 
strength 


Igrl = (28.1 ±1.4) MeV, 


(18) 


in perfect agreement with the number obtained in 
Ref. [13 from the decay /2 —tttt (compare also 
Ref. [57]), thus confirming SU(3) symmetry in this 
channel. 

The coupling of the 02 to pion and photon can be 
deduced from a Lagrangian (compare Refs. [2^129] 1 


Ctp^ = -^e^.„/3(T“^[/r,5^«A]), (19) 


where ± u’^F^'^u (omitting axial vector 

helds), = eQ{d^A'^ — d^'A^) is the electromag¬ 
netic field strength tensor, Q = diag(2/3, —1/3,...) the 
quark charge matrix, and we have neglected additional 
currents. Equation (HU) leads to the radiative decay 
width 


r{a2 Try) 


«-M^)5 

1607rF2 ^5^ 


( 20 ) 


which, compared to B{a 2 Try) = (2.68±0.31) x 10 
leads to 


\ct\ = (0.060 ± 0.004) GeV^ (21) 

If we combine the Lagrangians (fT51) and m with 
the tensor propagator iP^iv,pa{l)/[jn^^ — ^^), we can 
calculate the a 2 -exchange contribution Fa 2 (s,t,u) to 
77 —^ TT+TT^y. We hnd 


±■ 02 ( 3 , 1 , m) = g{s,t,u) + g{u,t,s), 
AeCTQT 1 


Q{sA,u) = 


V3F^ ml - s 


t — u + M/ — M 7 ; — 


2 {s + Ml){M^-Ml) 


( 22 ) 


which is completely fixed by experimental information 
up to an overall sign. 

A few remarks are in order concerning Eq. (I22E 
First, we can also perform an s-channel partial-wave 
expansion according to 

±(s,1,m) = '^Pl{zs)gi{s), 

I 


COS 9 s = Zs = 


s{t -u)- M2(M2 - Ml) 
(3-M2)A1/2(s,M2,M2)’ 


(23) 


which is the natural partial-wave expansion for jtt~ —> 
TT~g in terms of the scattering angle 0s- The partial- 
wave expansion of the s-channel 02 exchange amplitude 
Q{s.,t,u) then reads 


G{s,t,u) 

ffr(s) 

3zsgris) 


gr(s)±3rf(s), 

deCTffT ± ~ 

y/3Fl sml^ 


AecTgT 1 
\/3F3 ml^ - s 



s 

(24) 


Phrased differently, l/(s, 1, u) contains a nonresonant P- 
wave contribution (which has no 02 propagator) in ad¬ 
dition to the expected resonant D-wave. This is a well- 
known problem of higher-spin propagators; see e.g. the 
discussion in Ref. m- We cannot easily subtract the 
P-wave and use the P-wave alone, as Eq. (IM|) shows 
that both partial waves individually display an artifi¬ 
cial pole oc 1 /s, which is not present in the full ampli¬ 
tude (155)) . While a pole at s = 0 is not kinematically 
accessible in either of the two processes we consider in 
this article, it precludes a dispersive reconstruction of 
f-channel rescattering as discussed in the following sec¬ 
tion. We therefore retain the P-wave part in Eq. (15^ : 
its effect turns out to be numerically small. 

Second, we fix the sign of crgr in the following 
way. As pointed out in Ref. |4], the vector-meson con¬ 
tributions to 77 —>■ TT+TT^y determined in Ref. [3] can 
be rewritten, using the limit of a large number of col¬ 
ors (i.e., neglecting loop effects) and expanding the p 
propagators to leading order in the spirit of resonance 
saturation of chiral low-energy constants, as 


±(s, t, u) = Fri 


= F, 


rjTTTr'y 


i± 


i± 


3t 

2mj 


+ 0{m/) 


27772 j 


f2(t) + 0{m;^), (25) 


where we have used the approximation 

n 2 


n{t) 




t 


— t 




(26) 


In other words, Eq. (1551) predicts ~ l/(277i2) = 
0.84 GeV“^, a little more than half of the phenomeno¬ 
logical value aa « 1.52 GeV“^ (when using Eq. (fTIl 
for the definition oi a and not the pion vector form 
factor for a). We can now similarly expand Eq. (1551) 
to leading order in inverse powers of ml^. If we neglect 
the induced quark mass renormalization of the anomaly 
(proportional to M^, M^), we find the following esti¬ 
mate for the 02 contribution to an: 


Q-n — 


487r^CT^T 


±(0.46 ±0.04) GeV^ 


(27) 
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We shall see below that the true effect when including 
the 02 in a new extraction of the slope parameter an 
from data is significantly smaller, mainly due to curva¬ 
ture effects in the induced amplitude. Still, while effects 
in particular of excited p' resonances can be nonnegli- 
gible, we take the discrepancy between the p-induced 
slope and the experimentally determined value as an 
indication that the sign of the 02 contribution ought to 
be positive, 

CTflT = +|ct5t|- (28) 


The point of view of chiral perturbation theory allows 
us to further substantiate this choice. If we add the am¬ 
plitude (1^ . expanded to leading order in I/TOq^ 
with the sign as in Eq. (l28ll . as a further resonance 
saturation contribution to the one-loop representation 
of Ref. [3], the partial width E (?7 —^ 7r+7r“7) increases 
by about 7 eV, bringing the original prediction of 47 eV 
into even better agreement with the experimental num¬ 
ber (55 ±2) eV [35]. We will therefore work on from this 
hypothesis, and we give further hints below that data 
indeed suggests this to be the more likely solution. 

As a final remark, we will later insert a nonvanish¬ 
ing, energy-dependent width in the ai propagator in 
Eq. (1331) by hand. 


1 1 

- S - (s) ’ 

using the parametrization m 


(29) 


^02 ('®) — '^02 ^ ' Pi 

i=ri,p 


ma2 Qijs) T{qi{s)R) 
Vs QiimlJ T{q{mljR) ’ 


Qrj/p{s) — 


2 yi 


Tix) = 


9 -I- 3a;2 -I- x'^ ’ 


(30) 


which explicitly takes into account the 02 decays into 
final states Trrj and irp with relative branching frac¬ 
tions Prj = 0.17, Pp = 0.83, using the barrier factor 
R = 5.2GeV~^. The 02 is sufficiently far from the irp 
“threshold” that it seems a justifiable approximation to 
treat the p as a stable particle in this case. In contrast to 
using a constant width, this parametrization provides 
the correct threshold behavior of the imaginary part, 
as well as a reasonable phase above the resonance. 


2.3 Unitarization 

It is obvious that simply adding the tree-level 02 con¬ 
tribution (1331) to the original amplitude (ITUl) violates 
Watson’s theorem: we are missing the pion-pion rescat¬ 
tering on top of the a 2 -exchange graphs; see Fig.|3| The 



Fig. 3 a 2 contribution with pion-pion P-wave rescattering. 
The graph is supposed to comprise both possible charge as¬ 
signments inside the loop, corresponding to tree-level s- and 
M-channel 02 exchange; see Fig. [3] 


full dispersive solution that reinstates the correct phase 
relation in the t-channel P-wave is of the form 


R{s,t,u) = T{t) -I- g{s,t,u) + g{u,t,s), 

T{t) = f7(t)|A(l -I- ant) 

P 1'°° da; sin (5 (a;) ^ (a;) 1 
tt 74^2 a;2 |l7(a;)|(x - t) j ’ 

g{t) = ^J z^)[g{s,t,u) + g{u,t,s)]. 

(31) 


g{t) is the projection of the 02 exchange graphs onto 
the t-channel P-wave. It is given explicitly by 


0{t) 


_ Secrgr f _ 1 -|_1 


I "i; 


M^-t 


„ (m 2 +M 2 )(M 2 -M 2 ) 

+ 2t- 3M2 - - 11 


Q{y) = 


y = 


, 1 - y , y +1 
, y + —^ log —7 

o-t V 2 y-l 

2ml^ - M 2 - 2 M 2 + t 


Qiy), 


(32) 


g{t) contains a square-root singularity at t = 0, signal¬ 
ing the onset of the left-hand cut. As g{t) approaches 
a constant for large arguments t —>■ 00 , two subtrac¬ 
tions in Eq. m are sufficient, as the Omnes function 
behaves asymptotically as f7(t) ^ 1/t for S{t) —>■ tt. 
The number of subtractions therefore exactly reflects 
the original form in Eqs. m and cu¬ 
lt is easy to see that the full t-channel P-wave re¬ 
sulting from Eq. (ED), 

/i(t)=P(t) + e(t), (33) 


has the correct phase, while P(t) alone is subject to the 
inhomogeneous discontinuity relation 

discP(t) = 2i[p(t) -I- g{t)] sin(5(t)e“*'^l*ld(t — 

( 34 ) 
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The representation using the inhomogeneity Q{t) 
as input to the dispersive integral, preserves unitar- 
ity in the t-channel in the presence of left-hand cuts, 
which are approximated by resonance (here: 02 ) contri¬ 
butions. This is closely related to the methods used e.g. 
in Ref. [53] for 77 —>■ tttt, or in Ref. |3T] for semileptonic 
R-decays. We cannot easily apply an iterative proce¬ 
dure to determine left-hand cuts from right-hand cuts 
and vice versa, as done e.g. in the analysis of the closely 
related Primakoff process 77 r —>■ tttt as we do not 
have independent information on ttt] scattering phases 
at our disposal. 

Obviously, the 02 s- and w-channel exchanges will 
also generate nonvanishing projections onto F- and 
higher t-channel partial waves. These partial waves are 
real as long as we neglect pion-pion rescattering ef¬ 
fects in those higher waves, which is entirely justified 
for 7] —>• 7 r“'" 7 r “7 (and even for r]' —> 7 r“'' 7 r“ 7 ), given the 
smallness of the corresponding phases; compare the dis¬ 
cussion in Ref. m- However, even the real part of the 
T’-wave is entirely negligible: while in the chiral power 
counting, it is suppressed compared to the P-wave by 
another power of we have checked that kine- 

matical prefactors effectively suppress it by more than 
3 orders of magnitude in the physical decay region of 
rj —>■ TT+TT^y, and still by 2 for p' —>■ Tr+Tr^y. We will 
therefore discuss the comparison to decay data in the 
following section still in the approximation indicated in 
Eq. ©, using the P-wave only. 

3 Comparison to decay data 

3.1 p —>■ 7 r“'' 7 r “7 decay spectrum 

In this section, we compare the amplitude constructed 
in the previous section to the data on dP/dt as obtained 
by the KLOE Collaboration m- The decay distribu¬ 
tion was measured with arbitrary normalization, which 
has to be fixed independently from the branching frac¬ 
tion B{p —>■ 7 r+ 7 r“ 7 ) = (4.22 ± 0.08)%, as well as the 
total width of the p [26] . 

We first (re)fit the representation (fTUl) . (fT5]) . We ob¬ 
tain 

= (1.52 ±0.06) GeV"^ (35) 

where the error is only due to the statistical uncer¬ 
tainty in the data and neglects all the systematic ef¬ 
fects discussed in Ref. nni. The difference in the central 
value compared to a in Eq. (|3|) is due to employing the 
Omnes function instead of the pion vector form factor, 
as discussed above0 The quality of the fit is excellent, 

^In fact, if we construct the Omnes function from the phase 
of the pion vector form factor instead of from the tttt P-wave 



Fig. 4 Representation of the decay distribution rj 7r+7r“7 
from Ref. 1101 : see main text for details. The blue dashed 
curve shows the linear fit, while the full red curve includes 
the effects of 02 exchange in addition. The vertical dashed 
lines represent the limits of phase space at and M^. 

with a per degree of freedom of y^/ndof = 0.94. 
The subtraction constant A that, in this case, serves 
as an overall normalization of the amplitude, is H = 
(5.43 ± 0.12 =F 0.04) GeV“^, where the first error is due 
to the uncertainty in the integrated partial width, and 
the second due to the uncertainty in almost per¬ 
fectly anticorrelated with the latter. A thus seems well 
compatible with see Eq. (|T|). 

In Fig.m we plot the following observable, obtained 
from the data from Ref. unj: 

i.e. within the accuracy of the amplitude representa¬ 
tion without left-hand cuts, we expect to find P{t) = 
P(t)/P(0) = 1 ± ant. As the quality of the fit sug¬ 
gests, the linear curve (blue dashed) describes the data 
perfectly. 

Including the effects of 02 exchange (properly uni- 
tarized in the t-channel), the subtraction constant an 
in Eq. eiD has to be refitted to the data. We obtain 

= (1.42 ±0.06) GeV■^ (37) 

with x^/ndof = 0.90. The uncertainty of the 02 cou¬ 
pling constants induces an additional error in an of 
±0.01 GeV“^, which we will neglect in the following. 

phase shift 1241 as in Ref. 1331 , the central value of an reduces 
to 1.37GeV~^, rather close to Eq. ©• We disregard the ef¬ 
fects of varying the tt-k phase input in the following: they are 
compensated by corresponding shifts in oq to a very large 
extent, and they lead to insignificant uncertainties compared 
to other error sources. 
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The resulting fit is also shown in Fig. HI The reduc¬ 
tion in the value oi an compared to Eq. may seem 
surprisingly small, given the estimate of the 02 con¬ 
tribution to this parameter in Eq. (EZl). The reason 
is the curvature in P(t): in fact, the derivative P'{t) 
(which equals the constant an va. the simple fit) varies 
from P'{4.Ml) = 1.69 GeV^ to P'{M‘^) = 1.30 GeV^ 
within the decay phase space; outside phase space, we 
hnd e.g. .P'(l GeV^) = 0.46 GeV~^, and naive continu¬ 
ation to yet higher energies makes the derivative vanish 
and change sign around \/t = 1.25 GeV. It finally di¬ 
verges at f = 0 due to the square-root singularity. 


3.2 Impact on the 77 transition form factor 


As far as the phenomenological description of the rj —>■ 
7 r+ 7 r “7 decay data of Ref. m is concerned, the two 
amplitudes, with and without 02 effects included, are 
clearly equivalent: they describe the data equally well, 
and in fact, the two fit curves displayed in Fig. 0] de¬ 
viate from each other by less than 1 % in the whole 
decay region. This is different in the wider kinematic 
range of the similar decay rf —>■ 7 r“'' 7 r“ 7 , which we dis¬ 
cuss in Appendix A While the available data do not 


yet allow one to prefer one amplitude over the other 
in a statistically valid sense, the comparison of the ex¬ 
tracted subtraction constants an and an a'^ dehned in 
an analogous manner for rj' —^ 7 r“'" 7 r “7 seems to favor 
somewhat the decay amplitude including the curvature 
effects induced by the 02 . 

However, we have emphasized in the introduction 
that the decay amplitude 7 —> 7 r“'" 7 r “7 serves as a cru¬ 
cial input to a dispersive analysis of the 7 transition 
form factor [5], where the dispersion integral extends 
over a much larger range in energy (in principle, up to 
infinity). We therefore may expect to see a somewhat 
more significant deviation between the two amplitudes 
in there. 

We refer e.g. to Ref. [5] for all pertinent definitions 
concerning the singly virtual 7 transition form factor 
,0), which at small photon virtualities can be 
expanded according to 


F^^.t,(0,0) 


= 1 


Q" + OiQ^). (38) 


The slope parameter 6 ^ is divided into an isovector 
7=1 and an isoscalar 7 = 0 piece. The isoscalar part is 
small: employing uj + (p dominance together with data 
input on 777 yields ~ —0.022GeV“^ [5]. 

The slope is therefore almost entirely given by the 
isovector contribution, which in turn is dominated by 
7 r+ 7 r“ intermediate states; see Fig. [5] The correspond- 



Fig. 5 Two-pion cut contribution to the isovector part of the 
(singly virtual) t; transition form factor. Here, the gray circle 
denotes the t-channel rj 7r+7r“7 P-wave, while the white 
circle is the pion vector form factor. 


ing sum rule can be written as [S] 


da; 


5 (^= 1 ) =_^_ / 

J4M2 X 


atFris)h{x), 


(39) 


where F^(t) is the standard pion vector form factor, 
and we have written the dispersion integral with a cut¬ 
off instead of integrating to infinity. The 77 —>■ 77 
amplitude is obtained from the corresponding par¬ 
tial width by 


^77 = (40) 

Following Ref. [5], we vary the cutoff in the range 
_ {M^,... 2 GeV^}. With the decay amplitude (ITOl) . 
(fT^ . we find 

6^=1) = [{2.04 ... 2.22} ± 0.04„ 

± 0.02b ± O.Ol^y] GeV“^ (41) 

where the indicated range follows the range of cutoffs, 
and the errors are due to uncertainties in an^ the 
branching ratios for 77 —^ 7 r“'" 7 r “7 and 77 —^ 77 , and the 
pion vector form factor. For the latter, we employ the 
pion vector form factor parametrizations of Refs. PlIM] 
(or approximations thereof). Using, however, the par¬ 
tial wave /i (t) as in Eq. (1551) , the result reduces to 

= [{1.90...2.04}±0.04„ 

± 0.02b ± O.OIbv ± O.Ola^] GeV■^ (42) 

with the additional error due to the uncertainty in the 
02 coupling constants. That is, the slope is reduced by 
about 7%, a bit more than the combined error cited 
in Ref. [2, for a cutoff A^ fa IGeV^; this reduction 
is increased for higher cutoffs (due to the increasingly 
stronger curvature effects). A more detailed investiga¬ 
tion of 02 effects on the 77 (and 77 ') transition form 
factor(s), beyond the value of the slope at the origin, 
should still be pursued. 


^ While we have propagated the statistical error on from 
Eq. (1351) only in the rest of this paper, we here use the larger 
uncertainty ±0.13 due to systematic effects HOI in order to 
be consistent with the analysis in Ref. [5] . 
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Fig. 6 Radiative correction to 771 it rj due to photon- 
exchange diagram. 


4 Phenomenology for jtt —>■ tzt ] 

In the previous section, we have constructed an 77 —> 
7r+7r“7 decay amplitude including the leading left- 
hand-cut contribution, and have shown that this am¬ 
plitude describes the available decay data very well. As 
this representation includes the lightest resonance that 
can contribute in the nr] system, we are well equipped 
to now consider the crossed process 


7(fc)7r (P2) tt ipi)r]{q), 


( 43 ) 


which is described by the same amplitude as the decay 
process in Sect. [ 5 ] with Pi = —pi (using time-reversal in¬ 
variance). The Mandelstam variables are defined as be¬ 
fore, e.g. s = (pi+q)"^ denotes the total energy squared 
in the center-of-mass system, t = (p2 — Pi)^ is related to 
the pion momentum transfer etc. In particular, Eq. (ESI) 
is the natural partial-wave expansion in scattering kine¬ 
matics. 

The (polarization-averaged) differential cross sec¬ 
tion is given by 


da (s-M 2 )A 3 / 2 (s,M 2 ,M 2 ) 

df? 20487r2s2 


(1 - Zs)\J^{s,t,u)\^, 
( 44 ) 


from which one obtains for the total cross section 
(s-M2)A3/2(s,M2,M2) 


s) = 


10247rs2 


X J dzs{l- z'^)\J^{s,t,u)\'^ 

(s-M2)A3/2(s,M2,M2) 

7687rs2 

l5i(s)P + ^|52 (s)P + yl53(s)P 


( 45 ) 


where we have inserted the s-channel partial-wave ex¬ 
pansion (l 23 )) up to F-waves in the second step. 

As a cautionary side remark, we wish to point out 
that it has been emphasized in Ref. [ 35 ) for the similar 
process ■jn~ —>■ 7r“7r° that there is one significant effect 
due to radiative corrections, which is due to photon ex¬ 
change in the t-channel, compare Fig. [51 Translated to 



Fig. 7 Total cross section crfs) for 'yn~ n~ri. The blue 
band shows the cross section obtained from crossing the de¬ 
cay amplitude of Ref. [4]; the red band corresponds to the 
full amplitude including 02 effects. Finally, the yellow band 
displays the full cross section for the relative sign of the 02 
contribution flipped. The insert magnifies the near-threshold 
region. See main text on the error bands. 


the process under investigation here, the inclusion of 
this effect amounts to correcting the scattering ampli¬ 
tude in the form 


2e 

F(s,t,u) F(s,t,u) - y^77- 


( 46 ) 


Strictly speaking, the photon-exchange amplitude 
would have to be amended by form factor effects, in¬ 
cluding both the T] transition and the pion vector form 
factor; however, these corrections were shown to be very 
small in ^n~ —>■ n~n^ | 35 ] . The inclusion of the cor¬ 
rection (gSl) may be desirable if experimental data on 
77r“ —> n~ri become sufficiently precise in the future; 
we still neglect it for the following investigation. 

We show the total cross section in Fig. [T] We com¬ 
pare the cross section obtained from the decay ampli¬ 
tude in Ref. [ 1 ] by crossing to the full cross section in¬ 
cluding 02 effects. We find that dominance of t-channel 
dynamics holds roughly up to 7/s = 1 GeV, while above 
this value, the tensor resonance begins to dominate. We 
predict a peak cross section of about (12 ± 2)/rb, which 
is of a similar order of magnitude as the cross section 
of ‘jn~ —^ 7r“7r*^ at the p peak m- For completeness, 
we also display the cross section with the relative sign 
of the 02 contribution, see Eq. ESI, flipped (and all 
other parameters adjusted such as to best reproduce 
the 77 —>■ 7r“'"7r“7 decay data); we see that the transition 
from the near-threshold to the resonance region looks 
quite different, for reasons that will become transparent 
below. The uncertainty in the resonance peak is obvi¬ 
ously dominated by those in the 02 coupling constants 
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Fig. 8 First three partial waves for 777 “ —>• 7 r“r;. The moduli are shown in the normalization ii(Z + l)|gi(s)| for P-wave (blue 
bands), P-wave (red), and P-wave (green); bands with dashed borders refer to the analytic continuation of the decay amplitude 
in Ref. [3], while the full bands show the full result including 02 effects. The phase of the complete P-wave is represented by 
the red-striped band. All indicated bands combine the uncertainties in P (?7 ^ 7 r+ 7 r“ 7 ), an, and the 02 couplings ctQt- 


ct 5 Tj while near threshold, the errors coming from the 
total decay rate r[rj ^ 7 r+ 7 r“ 7 ) as well as an are more 
important. 

In the introduction, we pointed out that a naive con¬ 
tinuation of the 77 —>■ 7 r+ 7 r “7 decay amplitude Eqs. m 
and m would lead to a zero in the scattering ampli¬ 
tude 77 r“ — 7 > 7 r “?7 at t = —1/an. As s increases, this 
zero first appears in the differential cross section dF/dzs 
in backward direction, i.e. for Zg = —1. Given the form 
of the partial-wave expansion (ESI, 

F{s, t, u) = gi{s) + 3zsg2{s) + ..., (47) 

and assuming F- and higher partial waves are small, 
this will occur once the D-wave is one third the size 
of the P-wave, as long as relative phases are small. In 
our amplitude representation, the only imaginary part 
stems from the energy-dependent width of s-channel 02 
exchange; the P-wave phase is neglected, and all partial 
waves induced by t-channel exchange are obviously real. 
For better comparison and due to 

+ 1 )^ ( 48 ) 

we display the hrst three partial waves multiplied with 
5 Z(/ + 1) in Fig. [5] the intersection of P- and P-wave 
curves then gives an indication at the energy at which 
an additional zero in the angular distribution will oc¬ 
cur, with the precise position slightly modified by the 
small, but nonnegligible P-wave. We compare the full 
amplitude including the 02 to the continuation of the 
decay amplitude from Ref. [1] . The decisive observation 
is that including the 02 , the P-wave becomes more im¬ 
portant than the P-wave at even lower energies, around 


y/s = 0.9 GeV, where the phase is still tiny—we there¬ 
fore indeed expect to observe an almost perfect vanish¬ 
ing of the amplitude. To demonstrate that this is not 
trivially so, Fig. also shows what would happen with 
the opposite sign of the 02 contribution: negative inter¬ 
ference of s-channel 02 and t-channel exchange leads to 
a near-vanishing of the P-wave around 1.1 GeV (which 
is the cause for the rapid phase variation at that en¬ 
ergy) , and its rise toward the 02 peak only overtakes the 
P-wave once the phase is significant. As a consequence, 
no near-complete cancelation ever occurs at any energy. 

We wish to re-emphasize that there is no fixed rela¬ 
tion between the phase of our s-channel partial waves 
to m] scattering phase shifts according to a hnal-state 
theorem. As the corresponding ttt] phases are not the¬ 
oretically determined in the way the tttt [2^13611^ or 
ttK |38] phases are, unitarization using model phases 
seems to offer no significant improvement. Furthermore, 
the 02 is a largely inelastic resonance with respect to 
ng scattering anyway, with the dominant decay channel 
being irp [see Fq. dsni)], such that no simple version of 
Watson’s theorem applies, and any unitarization would 
have to implement a coupled-channel formalism. 

For illustration, we also show the resulting angular 
distribution at three sample energies y/s = 0.9 GeV, 
1.0 GeV, and 1.1 GeV in Fig. [51 indicating the transition 
between the threshold [P-wave dominance, da/dzg oc 
(1 — zf)] and the resonance region [P-wave dominance, 
dcr/dzs oc {1 — Zg)Zg]. The very different features of the 
different signs of the 02 are clearly visible. 

An observable that allows one to capture the key 
features of the effect discussed even in a comparably 
low-statistics measurement is to characterize the be- 
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Fig. 9 Differential cross sections Aa/dza for the three energies y/s = 0.9 GeV, 1.0 GeV, and 1.1 GeV (from left to right). The 
blue bands denote analytic continuation of the amplitude of Ref. [3], the red bands are our full predictions including a 2 effects, 
while the yellow bands show the same with the opposite relative sign of the a 2 contributions. The pronounced minima, very 
close to actual zeros, in backward direction in the red bands are clearly seen. 



Fig. 10 Forward-backward asymmetry according to 
Eq. (149 II . The color code is as in Fig. (9] 


havior in terms of a forward-backward asymmetry, 


_ Re ( 32 ) (gi + 33 ) , . 

“ 49,41 12 I 8„2 ’ 

9^1 + 5 152 I +733 

where in the second line we have neglected all partial 
waves beyond F-waves, as well as imaginary parts of 
the P- and the F-wave. As shown in Fig. [Tni both the 
continuation of the decay amplitude without 02 [ 1 ] and 
our full model with the preferred relative sign for the 02 
display a very large positive asymmetry, peaked just be¬ 
low y/s = 1 GeV for the full model; for the opposite 02 
sign, the asymmetry is small near threshold, and sub¬ 
sequently even turns negative. An experimental verifi¬ 
cation of this asymmetry would therefore confirm that 
our description of the decay amplitude including the 
02 , and the resulting consequences for the ry transition 
form factor, are indeed reasonable. 


5 Summary 

In this article, we have studied the effects of the 02 
tensor meson on the decay 77 —>■ Tr+Tr^y as well as the 
analytic continuation of the decay amplitude for the 
scattering process 771 “ —>■ Tr~r]. We have included the 
D-wave nr] resonance as a left-hand cut structure of a 
dispersive representation that obeys the correct final- 
state phase relation for the tt+tt” P-wave. While the 
decay spectra measured by the KLOE Collaboration 
can be described equally well with and without the 02 
effects, there seems to be an indication for better con¬ 
sistency of the subtractions constants when comparing 
to the similar decay rj' —>■ Tr+Tr^y. The slope parame¬ 
ter of the resulting ry transition form factor is reduced 
by about 7% in the dispersive integral up to 1 GeV^ 
compared to a previous analysis [5] . 

We have predicted different observables for the ry 
production reaction 'yn~ —>■ n~r] at energies up to the 
02 resonance. The peak cross section is predicted to be 
( 12 ± 2 )/ib, similar in size to the ■jn~ —>■ n~n^ cross sec¬ 
tion in the p peak m- Fixing the relative sign of the 
02 to the more likely solution from decay phenomenol¬ 
ogy, we hnd an interesting P-P-wave interference ef¬ 
fect, leading to almost perfect zeros in the differen¬ 
tial cross section, and a very strong forward-backward 
asymmetry in the energy region between threshold and 
the 02 peak. These predictions provide strong motiva¬ 
tion to study the corresponding Primakoff reaction e.g. 
at COMPASS, which may help to further scrutinize the 
physics of light mesons relevant for hadronic corrections 
to the muon’s anomalous magnetic moment. 
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Appendix A: rj' —)• 7r7r7 


The formalism of dispersively analyzing decay data 
with a final-state P-wave pion pair was also applied 
to 77 ' —?► 7 r+ 7 r “7 in Ref. [3], analyzing data by the 
Crystal Barrel Collaboration [52] • We do not intend to 
make a prediction for the crossed Primakoff reaction 
77 r — nr ]'—the threshold is too high, too close to the 
02 resonance tail to still find traces of the t-channel 
exchange, and the number of inelastic (subthreshold) 
channels probably too large to be ignored. However, we 
expect the impact of left-hand cuts in the decay process 
to be much stronger over the wider kinematic range ac¬ 
cessible in the rj' decay, i.e., the curvature effects that 
are rather moderate in the 77 decay in Fig. |3| should be 
much more visible in that case. Furthermore, the decay 
77 ' —>■ TT+TT^y is about to be remeasured with increased 
precision by BESIII (see [30] for spectra not yet cor¬ 
rected for acceptance), such that a prediction for the 
curvature of the spectrum (after dividing out the uni¬ 
versal Omnes factor) ought to be very timely. 

To determine the 02 contribution to 77 ' —>■ 7 r+ 7 r“ 7 , 
we first need to fix the 02 —> nr]' coupling constant. We 
can first do this by naively defining a coupling g'j, with¬ 
out a Lagrangian, just through the analogous relation 
to Eq. (fT71) 


P(a 2 nr]') 


{gr? 

1807rP4 rnl^ 


(A.l) 


which, with B{a 2 —>■ nr]') = (5.3 ± 0.9) x 10 yields 


Iffrl = (25.5 ± 2.3) MeV. 


(A.2) 


If we attempt to explain this value based on the sin¬ 
gle Lagrangian term CSI), we first need to amend the 
pseudoscalar field </> by the SU(3) singlet 770 , 

TfO + V2n+ 

V2n- -n° + 



^+.... (A. 3 ) 



Fig. 11 Representation of the decay distribution rj' 
7 r+ 7 r “7 from Ref. ESI; see main text for details. The blue 
curve shows the linear fit, including the gray band for the fit 
uncertainty. The red curve with the yellow band includes the 
effects of a 2 exchange in addition. The vertical dashed lines 
represent the limits of phase space at 4M.^ and . 


that we are confident the difference can be explained 
by higher-order terms. In particular, we can safely con¬ 
clude that the sign of g'rp in Eq. (IA.2I1 agrees with the 
one of (7 tI1 

In Fig. [Til we display the observable P"' {t) defined 
in strict analogy to Eq. (1551) . comparing the data of 
Ref. [39] to fits with a linear parametrization, as well 
as including effects of the 02 . Due to the rather large 
error bars, we show the ht results as bands, not just the 
best fit. The linear fit leads to a slope parameter 

a'jj = ( 0.6 ± 0 . 2 ) GeV■^ (A. 5 ) 

with a reduced of x^/ndof = 1.23. It was argued 
in Ref. [3] that in the limit of a large number of col¬ 
ors, an = a'^ should be expected, so the slopes of the 
polynomial would agree for 77 and 77 ' decay. Compar¬ 
ing Eqs. ill) and ESI), phenomenology seems rather 
at odds with this prediction. However, including the ef¬ 
fects of the 02 in the amplitude representation, we find 
much stronger curvature effects than for the 77 decay as 
anticipated, with the residual slope fitted to be 


In Eq. (1151) . we have simply identified the 77 with the 
octet field 773 ; if now we assume a simple, single-angle 
7777 ' mixing scheme, 

I77) = cos 6I|778) - sin6»|77o), 

I77') = sin6»|778)-b cos6»|77o), (A. 4 ) 

we can explain the ratio of the couplings by a mixing 
angle oi 9 = (—12.4±2.7)°, which is somewhat smaller 
than the standard value 9 ~ —20°, but close enough 


a'^ = (1.4 ± 0.4) GeV“^ (A. 6 ) 

now with x^/iidof = 1.38. The additional uncertainty 
due to the 02 couplings is ±0.1GeV“^. In this case, 
the fit quality becomes slightly worse (overall better 

®As a side remark, we point out that fixing an effective 
02 ^ np coupling constant from the known branching frac¬ 
tion 6(02 np) should also allow us to include 02 effects in 
the decays rj' —>• 47 r I41II42I , thus going further beyond vector- 
meson dominance. 
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fits are essentially precluded by the third-to-last data 
point at = [800, 825] MeV); however, u'q is now in 
markedly better agreement with the value found for a a 
in Eq. (EH). P'^ (t) can be approximated in the decay 
region 4M^ < t < M^, by a quadratic polynomial 

+ + (A.7) 

to about 1% accuracy. The 02 contribution predicts the 
curvature to be = (—1.0 ± 0.1)GeV“^. As a side 
remark, we can also take this result as another strong 
indication on the correctness of the sign of cxgx and 
CTg'rp'- a negative sign would lead to a residual slope o!q 
of (0.06 ± 0.12) GeV^ 

A more rigorous test of the decay spectrum pre¬ 
dicted here, with higher-statistics data from BESIII, 
would be extremely welcome. 
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